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We study dynamical behaviors of one-dimensional stochastic lattice gases with repulsive interac- 
tions whose span can be arbitrary large. We endow the system with a zero-temperature dynamics, 
so that the hops to the empty sites which would have led to the increase of energy are forbidden. We 
assume that the strength of interactions sufficiently quickly decreases with the separation between 
the particles, so that interactions can be treated in a lexicographic order. For such repulsion pro- 
cesses with symmetric nearest-neighbor hopping we analytically determine the density-dependent 
diffusion coefficient. We also compute the variance of the displacement of a tagged particle. 
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I. INTRODUCTION 

Lattice models play a central role in equilibrium and 
non-equilibrium statistical mechanics [T]-[5]. Many lat- 
tice models admit the lattice gas representation in which 
each lattice site is either occupied by a particle or empty, 
and particles in nearby sites interact according to some 
specified laws. In equilibrium statistical mechanics, the 
goal is to compute a free energy as a function of den- 
sity and temperature, and to understand the properties 
of phase transitions which often occur in lattice gases in 
d > 2 dimensions. In non-equilibrium setting, namely for 
lattice gases supplemented with stochastic dynamics, in- 
teresting dynamical behaviors and rich non-equilibrium 
steady states arise already in one dimension [3HS]- 

Here we consider one-dimensional lattice gases which 
are endowed with stochastic dynamics. We analyze var- 
ious stochastic zero-temperature dynamics which forbid 
energy-raising moves. Many our results are exact in the 
finite setting, viz. for N particles on a ring with L sites, 
but we focus on the thermodynamic limit, N — > oo and 
L — > oo with p = N/L kept finite, which is particularly 
important. We consider lattice gases with exclusion prop- 
erty (each site contains at most one particle) , so that the 
density lies in the range < p < 1. 

Among lattice g simple exclusion process (SEP) 

is the most tractable. In the SEP, particles undergo 
a nearest-neighbor hopping and thereby interact only 
through the exclusion property (no more than one parti- 
cle per site) . A rich set of behaviors has been discovered 
for the SEP, see [5HH] and references therein. Here we 
investigate a repulsion process (RP) where in addition to 
exclusion, there is an energy cost when particles occupy 
adjacent sites. We postulate the following dynamics: 

1. If an attempted hop (to a neighboring empty site) 
would not increase the number of nearest-neighbor 
pairs of particles, it is always performed. 

2. If a hop would increase of the number of nearest- 
neighbor pairs of particles, it is never performed. 



One can associate a Hamiltonian 

L 

"K = n mi+i (1) 
i=i 

with the repulsion process. Here = 1 if site i is oc- 
cupied, and rii = if it is empty. The above dynamics 
is based on the exchange, rii o rij+i. This exchange 
cannot occur if it would lead to the increase of the num- 
ber of nearest-neighbor pairs of particles, and thereby to 
the increase of energy. At zero temperature, energy rais- 
ing moves are forbidden, so the dynamics underlying the 
repulsion process is a zero-temperature dynamics associ- 
ated with Hamiltonian ([TJ. 

The RP is a very basic lattice gas model and it nat- 
urally arises in various contexts. For instance, the KLS 
model proposed by Katz, Lebowitz, and Spohn |10j turns 
into the RP at a certain point of the parameter space. 
The KLS model has played a major role in elucidat- 
ing key conceptual issues in non-equilibrium dynamics 
[TTMl4| . The one-parameter version of the KLS model 
can be interpreted as the repulsion process with finite- 
temperature dynamics, and it reduces to the RP at zero 
temperature and to the SEP at infinite temperature. A 
number of our results for the RP admit a generalization 
to the KLS model [TH [T5] , but the tools are very different 
and one can go much further for the simpler RP model. 

The PR also arises in the context of the Ising model on 
the square lattice supplemented with zero-temperature 
spin- flip dynamics |16j . (This was the original motivation 
for this study.) In a study of the evolution of a quadrant 
of minority phase, one can map the interface dynamics 
onto a one-dimensional lattice gas. In the simplest case 
of nearest-neighbor interactions, the lattice gas is merely 
the SEP (with symmetric hopping), adding next- nearest- 
neighbor interactions leads to the RP, adding interactions 
with 3 rd nearest neighbors leads to a repulsion process 
with next-nearest-neighbors, viz. with zero-temperature 
dynamics associated with Hamiltonian 

L L 

Jf 2 = J\ ^njn.i + i + J 2 ^nin i+2 (2) 
i=i i=i 

In this paper, we first study the RP with Hamiltonian 
|l]). We then show how to extend the results to the RP 
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with Hamiltonian ([2| with arbitrary J\ > J2 > 0. Fi- 
nally, we outline the extension to an infinite family of 
lattice gases class, namely to a class of generalized RPs 
with interactions of arbitrary finite range. 

The rest of this paper is organized as follows. In the fol- 
lowing section we present some of our main results with- 
out derivations. In Sect. |III[ we determine the structure 
of steady states and compute correlation functions for 
the RP with Hamiltonian ([I]). In Sect, 
various correlation functions. The resu 



IV 



we compute 



ts of Sects. HIH - 
|TV| are then applied to the determination of the diffu- 
sion coefficient and to probing the self-diffusion process 
(Sect. [V|. In Sect. VI we investigate the RP with Hamil- 
tonian (121). The influence of repulsive interactions with 
longer span is analyzed in Sect. |VlH We conclude with a 
discussion (Sect. VIII). Some details of the calculations 



are presented in Appendices. 



II. MAIN RESULTS 

We are chiefly interested in the RP with symmetric 
hopping. Our methods equally apply to the driven case 
(when the hopping is biased) since for the RP on the 
ring the structure of the steady states is independent on 
the bias and its strength. Hence whenever we present 
results for the driven RP, we limit ourselves to the totally 
asymmetric RP in which only hopping to the right are 
allowed (the rate of allowed hops is set to unity as in the 
symmetric version). 

The structure of steady states is remarkably simple: 
They are characterized by the maximal number of is- 
lands. When p < 1/2, each occupied site is delimited at 
both ends by vacant sites. Thus after a transient period 
the system will be in a configuration like 



• oo»o»o»o»ooo»o»o»oo»o 



(3) 



Hereinafter • denotes a particle and o denotes a vacancy. 
(In (§ we have illustrated the maximal-island configu- 
ration on the ring with N = 9 particles and V = 13 
vacancies.) The system will then forever wander on the 
phase space of such maximal-island configurations. This 
phase is clearly connected as follows e.g. by noting that 
each maximal-island configuration can evolve into the 
maximal-island configuration containing the longest pos- 
sible string of alternating particles and vacancies, namely 
• o • o • • • • o, followed by the string of vacancies. 

Similarly in the high-density region, p > 1/2, in every 
steady state configuration each vacant site is delimited 
at both ends by occupied sites. Low- and high-density 
regions are treated in the same way and the results are 
related via the mirror symmetry, pf> 1 — p. Hence it suf- 
fices to analyze one of the two regions. In the following, 
we always consider the low-density region if not stated 
otherwise. 

The remarkable feature of maximal-island configura- 
tions is that they all occur with equal probabilities. Once 
one establishes these assertions about the structure of 



steady states configurations and their probabilities, the 
task of computing the average current, the correlation 
functions, etc. essentially reduces to a rather straightfor- 
ward combinatorial analysis. The emerging results are 
exact for finite systems, so the thermodynamic limit is 
taken at the end and there are no preliminary assump- 
tions which are often (tacitly) implied in the analysis of 
infinite-particle systems. 

Here are some of our findings. The probability of each 
maximal- island configuration is 6 , where 6 is the total 
number of maximal-island configurations with N parti- 
cles and V vacancies that can be arranged on a ring of 
size L = N + V: 



e = 



V - 1 
N - 1 



L (V — 1 
N \N - 1 



(4) 



in the low-density region N < V. For the totally asym- 
metric RP, the steady state current J is 



J 



fV-2\ 
KN-1J 



N V -N 



Q + (n-\) lv ~ 1 

In the thermodynamic limit 



J(P) 



(irfcl) I < p <i 



(5) 



(6) 



where the mirror symmetry allows to deduce the current 
in the high-density region. 

The pair correlation function is given by 

L<72J ft A (V-t+q-l\ 

{mni+e) =2^ 1 7v\ TF^TV ( 7 ) 



q=l 



9 ~ 1 / (n) + (n-1) 



In the thermodynamic limit the pair correlation function, 
more precisely the connected pair correlation function, 
has a strikingly simple form 



{niTij) c = (riirij) - p 2 = p(l - p) 



1-P 



valid for all i and j. Similarly 



(niiijjc = p(l-p)[- 



P 



\3~i\ 



(8) 



(9) 



in the high-density region, p > h. Thus the connected 
pair correlation function exhibits a pure exponential de- 
cay modulated by an oscillating sign. 

We have also computed higher-order correlation func- 
tions. They exhibit remarkable factorization properties, 
e.g. in the low-density regime the three particle correla- 
tion function reads 



(ninjUk) 



(ninj)(njn k ) 



(10) 



for all i < j < k. 
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Lattice gases with symmetric hopping usually exhibit a 
universal hydrodynamic behavior, namely on large spatial 
and temporal scales they are governed by a non-linear 
diffusion equation [3] 



dp d 
dt dx 



D( P ) 



dp 

dx 



(11) 



The diffusion coefficient D(p) represents the spread of 
disturbances on an equilibrium background with uniform 
density p. The diffusion coefficient depends on the density 
in most lattice gases thereby making (111 a non- linear 
partial differential equation. Thus all the microscopic 
details of lattice gas dynamics are absorbed into a single 
number, the diffusion coefficient D(p). 

The diffusion description has been justified for numer- 
ous lattice gas models [3J 0] , but the diffusion coefficient 
D(p) has been computed in a very few models, mostly in 
the situations when the diffusion coefficient is indepen- 
dent on the density. We have found that for the RP 




< p < \ 
\<p<\ 



(12) 



The limiting behaviors, D(p = 0) = D(p = 1) = 1 
in (12), can be understood without calculations. When 



p — > 0, it suffices to analyze a single particle in the empty 
lattice. This particle diffuses with D = 1 (recall that we 
have set the hopping rates to unity) . Similarly when the 
density is very close to the maximal density p = 1, we 
can consider a single vacancy in otherwise fully packed 
lattice. To establish D(^) — 4, let us slightly perturb 
the p — \ state. There is just one half-filled equilibrium 



configuration, . . 
perturbation 



• o»o»o»o. 



Consider its smallest 



•o»o»o»»o»o»o»o 



with a single doublet. In this configuration both particles 
forming the doublet can hop. A single hopping event, say 
to the right, leads to 



o • • o • o 



o • o • • o 



Therefore the doublet effectively undergoes a symmet- 
ric hopping on distance two, so its diffusion coefficient 
is 4 times larger than the bare diffusion coefficient of a 
particle. 

Apart from diffusion, we analyze self-diffusion: We tag 
a particle and probe its long time behavior. For the one- 
dimensional RP (and other lattice gases with nearest- 
neighbor hopping and exclusion), the mean-square dis- 
placement of the tagged particle grows as 

(X\t))=D(p)Vt (13) 

We have found that the amplitude D(p) is given by 



Vtt p 2 



Ul-2p) 
I (1 - p)(2p ■ 



< p < \ 
\<p<\ 



(14) 



for the RP. In contrast to the diffusion coefficient which is 
maximal at half-filling, the self-diffusion amplitude van- 
ishes at half-filling. 

We now turn to the generalized repulsion processes 
(GRPs) which are characterized by the Hamiltonian 



— j\ 



UiTli+l 



i=l 



(15) 



The coupling constants J\ , . . . , J m are assumed to be pos- 
itive, so that interactions between particles separated by 
distance < m are all repulsive; particles separated by 
larger distance do not interact. In addition we assume 
that the coupling constants obey the set of constraints 



Jk > Jk+l 



Jrr 



k = 1, 



1 



(16) 



An attempted hop to the neighboring empty site is again 
performed if the energy is not raised; otherwise, it is re- 
jected. The constraints (16) simplify the analysis since 



they allow us to treat interactions in a lexicographic or- 
der. Namely, if the attempted hop would lead to the 
decrease of the number of nearest-neighbor pairs, it is 
always performed; if the number would increase, it is 
never perform. If the number would remain the same, one 
should check how the number of next-nearest-neighbors 
would be affected. If the attempted hop would lead to 
the increase or decrease of this number, the fate of the 
hop is settled (it is rejected or accepted, respectively); 
otherwise we go to the third level. 

There are m relevant low-density regions with different 
behaviors 



< p < 



m+l ' 



<p< 



\<P<1 



(17) 



In each region, a certain class of maximal-island config- 
urations constitutes the phase space of steady states; in 
each class, the maximal-island configurations occur with 
equal probabilities. The corresponding maximal-island 
configurations are specified by the allowed range for the 
length s of islands of vacant sites: 



< p < 
l 

m+ 



T < P< 



<p< 



s > m 

s = (m — 1, m) 
s = (m — 2, m — 1) 



(18) 



etc. 



ending with region | 
and 



< 



where s 
l 



(1,2). 
are new, other 



< P - 2 

Regions < p < ^ and ^ < p < 
regions already appear in the analysis of the model with 
repulsion range m — 1. Thus in the \ < P < \ density 
range, the behavior stabilizes when m > k. 

The total number of admissible maximal-island con- 
figurations is a generalization of Q in the region with 
smallest density 



e = 



L fV 

N 



(m-l)N- 1 
N — 1 



0<P<^TT 



(19) 
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In the remaining low-density regions 

L ( N \ 1 

N \V - (k- l)N 



e = 



k+i <p< k 



(20) 



where k = 2, 3, . . . , m. 

Most of previous results for the RP (which corresponds 
to m = 1) admit an extension to the GRPs. Formulas 
for the single low-temperature region must be replaced 
by their analogs for the m low-temperature regions (171, 



but extensions tend to be straightforward. For instance, 
for the totally asymmetric GRP, the steady state current 
in the low-density regions is given by 



p[l-Qn+l)p] 
X—mp 

[(fc+i)p-i][i-fc P ] 



0<p< 

i 

fe+i 



1 <p<1 



(21) 



where k = 2, 3, . . . , m. 

The behavior of the diffusion coefficient characterizing 
the GRP with symmetric hopping is even simpler, for 
every m > 2 there are just 4 different regions: 



(1 — mp)~ 

„-2 



D(p) = 




^TT <P< 3 

- < o < m 

m+1 " 



(22) 



The self-diffusion amplitude in the low-density region is 



p 2r D = 



p[l - (m + l)p] 
[(fc+l)p-l][l-fcp] fe+1 



0<P<^TT 
1 <P<1 



where again k = 2, 3, . . . , m. The behavior in the high- 
density regime is obtained by replacing p — > 1 — p on the 
right-hand side. For instance, 



D 



_2_ [fc- (fc + l)p][fcp- (fc- 1)] 
sJtx p 2 

k 



when < p < j^xj, where fc = 2, 3, . . . , m. 

The results are particularly neat for the GRPs with in- 
finite interaction range. The constraints ( 16 ) are obeyed 



when e.g. the coupling constants decay exponentially 
with distance, Jk — 8 k with 8 < |. For such GRPs with 
totally asymmetric hopping, the current is given by 



J(p) 



[(fc+i)p-i][i-fc P ] 

[k-(k+i)p][k P -(k-i)] 
i-p 



kh <p< \ 

h T k <P<Wi 



(23) 



with k = 2, 3, . . .. For the same model with symmetric 
hopping, the diffusion coefficient is 



D(p) 



pY 



< p < I 

\ <p< 1 



while the self-diffusion amplitude is given by 



D = 



[(fc+i) P -i][i-fcp] 
p 2 

[fc-(fc+i)p][fcp-(fc-i)] 
p 2 



kTi <P< S 
' 5 ¥ 1 <P< ' 



(24) 



(25) 



with fc = 2, 3, Comparing ( 24 ) with ( 12 ) we see a sur- 



prising formal duality between the density-dependence 
of the diffusion coefficient for the RP with infinite-range 
interaction range and for the simplest RP with nearest- 
neighbor interactions. 



III. STEADY STATES 

In this section we explain the structure of steady states 
outlined in Sect.[llJ and we derive the announced results 
for the number of maximal-island configurations, Eq. Q, 
and for the steady state current in the totally asymmetric 
case, Eq. (5). 

We put our stochastic lattice gas on the ring and let 
it evolve. The initial condition will be eventually forgot- 
ten. To appreciate the emergence of configurations like 
([3| we notice that in each hopping event the number of 
islands either increases or remains the same; the decrease 
of the number of islands would imply the increase of the 
energy which is forbidden. Thus the system eventually 
reaches a state where the number of islands is maximal 
and then it wanders on the phase space of such maximal- 
island configurations. To prove that all maximal-island 
configurations are equally probable we denote by P{C) 
the probability of being in maximal-island configuration 
C and write the stationarity condition 

P{C) T (C -^ C ") = J2 p {C")T(C" -> C) (26) 

c C" 

where T(C — > C) is the transfer rate from C to C . 
Since T = 1 if the evolution is allowed and otherwise, 
we need to count the number of ways into and out of 
a configuration. Consider for concreteness the totally 
asymmetric version of the RP. The basic hopping event 
o • oo — > o o »o shows that the evolution out of C occurs 
at the left edges of islands of vacancies of length > 2. 
The same basic process shows that the evolution into 
C is counted by the number of right edges of islands of 
vacancies of length > 2. Therefore J2c> T ( C ~> C ') = 
J2 C » T ( C " -> C ) implying that if P(C) are equal for all 
configurations, Eq. (261 is satisfied. 



To determine C we proceed in the same way as in |30j . 
We pick up an arbitrary site as a first site and represent 
configurations as strings. Configurations in which the 
first site is empty look like 

N particles 



V— 1 vacancies 



with N particles distributed among V possible positions 
(each such position is denoted by 4-)- The number of 
such configurations is (J\ . Similarly configurations with 
occupied first site look like 



N—l particles 



• o 4,o|o|o4,o4,o|o|,o|, 



fc+1 



V— 2 vacancies 



5 



and the number of such configurations is ( N _ 1 )- Com- 
bining these two contributions we arrive at Q. 

To compute the steady state current J we count the 
number of configurations of the type 



AT— 1 particles 



o • o o loioio^oloioj, 



V — 3 vacancies 



and divide it by the total number of maximal-island con- 
figurations given by Eq. Q. This leads to Eq. 



IV. CORRELATION FUNCTIONS 

In this section we compute correlation functions. The 
results for the pair correlation functions, Eqs. ([7|-([9|, 
have been announced in Sect. |nj We also establish ( flp| ) 
and show how to compute other higher-order correlation 
functions. 

Let us first compute the pair correlation function 
(niUj) in the low-density regime. In the steady state, the 
pair correlation depends only on the separation £ = j — i 
between the sites. Obviously 



(n?) = (rij) = p, (riin i+1 ) = 



(27) 



The pair correlation function becomes non-trivial when 
1 = 2. This describes two particles separated by one site 
(which is necessarily empty in the low-density regime). 
We need to count the maximal-island configurations of 
the type 



N— 2 particles 



(28) 



V— 3 vacancies 



One gets ( N _ 2 ) , from which 



HiH+2/ 



\N-2) 



O + ( v N -\) 



1 



(29) 



Generally to compute (riini+e) let us consider configu- 
rations with q — 1 particles and £ — q vacancies between 
sites i and i + £: 



q—l particles 



N—q—1 particles 
A 



o»ojo|ojojoto|ojo|o|o|ojoj (30) 



-q vacancies 



V —i+q—2 vacancies 



Particles at sites i and i + £ are explicitly presented in 
(30), possible locations of other particles are shown by |. 



The number of configurations of type ( 30 ) is 

-q-l\ fV-l + q-X* 
q-1 )[ N-q-1 



(31) 



The first binomial factor accounts for placing q—l par- 
ticles into £ — q — 1 possible positions inside the string 



+ 1, . . . ,i + £). The second binomial factor describes 
all possible placings of N— q—l particles into V — £+q — 1 
admissible positions outside the string (i,i + l, . . . ,i + £). 
Taking into account that q varies from q = 1 to q = 
we arrive at (17]). In the thermodynamic limit, Eq. ^ 
becomes 



{n t n i+e ) 



E 

9=1 



£-q-l 

q-i 



P 9+1 (l-2p) 



l-2q 



(32) 



In particular, for 3 < 
reads 



(i - p y-* 

< 7 the pair correlation function 



(nin i+3 ) 


1 - 


2p 




P 1 


(1- 


P? 






(1- 


2pf 




P 1 


(1- 


-pf 


(1-P) 2 


(nin i+5 ) 


(1- 


2pf 


2p(l - 2p) 


P 2 


(1- 




(l-p)3 


(ntrii+e) 


(1- 


2pf 


3p(l - 2p) 2 


P 1 


(1- 


-pf 


(1-P) 4 


(n^i+r) 


(1- 


2pf 


4p(l - 2p) 3 


P 2 


(1- 


-pf 


(1-pf 



(33) 



(1-pf 
P 2 (l-2p) 
(1-P) 4 



The remarkable feature 
lation functions factorizc: 



of the SEP is that the corre- 
(niHj) = (m) (rii) = p 2 , etc. 



We have written ( 33 ) in a way that emphasizes that for 
the RP, the pair correlation function does not factorizc. 
Despite this lack of factorization, the connected pair cor- 
relation is given by a simple formula ([8]). This is easily 
verified for £ < 7 using explicit expressions (27), (29), 
and ( 33 1 . To derive Eq. j8l) from Eq. ( 32 ) in the general 



case, consider first the situation when the distance £ be- 
tween the sites is even. Writing £ = 2k + 2 and q = 1+p, 
we transform the sum on the right-hand side of ( 32 ) into 



? / 1 rt \ 2k k 

p 2 1- 2p s 



P(l - P) 



i-p v i-p j f^ \ p r ' (i-2p) 2 

The last sum can be computed for arbitrary k 



k 

E 

p=0 



2k -p 
V 



1 



1 



1 - 



V4z + 1 
1 



where 



X, 



V4.x + 1 
2x + 1 ± y / 4x+ 1 



Combining all these results we establish Eq. (8| in the 
situation when the distance £ — j — i between the sites is 
even. The situation when the distance between the sites 
is odd is treated in a similar manner. 

We now compute higher-order correlation functions. 
We begin with the three particle correlation function. 



(> 



We shall use the shorthand notation (123) = (n^n^n^). 
Let i = ii > and m = 13 — 12 > be the separations 
between the particles. Configurations 



9 -i 



N—q—r—l 



oioj,o|o|o«o|o|o|oto|o|o|.o|oj.o|oj 

i~q m—r V—l—m+q+r — 2 

represent arrangements with q — 1 particles between sites 
ii and 22 and r — 1 particles between sites Z2 and 13. The 
total number of such configurations is 

— q — 1\ fm — r — 1\ /V — £ — m + q + r — 1 
so that 



[£/2\ y m /2\ 

(123) EE 

q=l r=l 



-9-1 

9-1 



m — r — 1 
r- 1 



V iV-g-r- 1 / 

x O + (53) 

which in the thermodynamic limit becomes 

(123) = L ^ J (£-q-l\ ffl(l-2p)'-** 

p ^ v 9-1 y (1-p)^ 9 

Lm/2J 'm - r - 1\ p r (l - 2p) m - 2 '- 
r - 1 / (1 - p)" 1 -'' 



E 



The first sum on the right-hand side of the above equa- 
tion is equal to p~ 1 (l2) 1 while the second sum is equal 
to p^ 1 (23). (Here (12) = (n^n^) and (23) = {n i2 n i3 ); 
also p = (1) = (2) = (3).) Thus a seemingly cumbersome 
expression for the three particle correlation function can 
be re-written as 



(123) = 



(12) (23) 
(2) 



(34) 



which is Eq. ( |To| announced in Sect. [IT] Equation (34 1 is 
reminiscent to the Kirkwood's superposition approxima- 
tion [T7] which is often used in liquid theory [TH1 [H] . In 
the present case, however, it is an equality rather than 
an uncontrolled approximation. 

For the RP, the correlation functions do not factor- 



( 34 1 the rij a factors which appear twice in the numerator 



would be balanced by appearance in the denominator 
thereby agreeing with a single appearance on the left- 
hand side of (135). 



By inserting (8b into (35) we recast the higher-order 



correlation functions into a sum of exponential factors. 
The simple form of Eq. ^ suggests that the connected 
correlation functions defined via [18 



M a ) c = (Kj - /£>)••• (n ia - p)) 



(36) 



may have a simpler form. The three-particle connected 
correlation function 

(123) c = (123) - p [(12) + (23) + (13)] + 2p 3 



simplifies to a single exponential 



(123)o = p(l-p)(l-2p) 



1-P 



(37) 



In Eq. (37) we assume that ii < i 2 < 13; when some 



indexes are equal, the predictions remain correct. Similar 
remarks apply to following formulae with more indexes, 



e.g., to Eqs. (38), (All, (A2) 



The four-particle connected correlation function re- 
duces to a combination of two exponential terms 



(1234) c = (l-2p) 2 p(l-p) 



P 



+ W-p)?[- 



1-p 



1-p, 



(38) 



Equations (|37|)-(|38j indicate that the number of ex- 
Explicit formulae 
in Ap- 



ponential terms increases with order 
for (12345) c and (123456) c [Eqs. dATJ and (TA2 



pendix |A] confirm this assertion. Generally tne num- 
ber of exponential terms appearing in the p— particle 
connected correlation function is equal to the Fibonacci 
number F p -i (Appendix |A|. Hence the number of expo- 



1] 



1.618 



nential factors varies as g p , where g = | ( v5 - 
is the golden ratio. In comparison with 2 P ~ 1 exponen- 
tial factors which appear in ordinary correlation function 
(12 . . .p), there is a huge cancellation for the connected 
correlation function (12 . . .p) c , although the growth with 
p remains exponential. 

The behavior in the high-density region is established 
by replacing n, — > 1 — n,; and p — > 1 — p. For instance, 



ize as we already know. However, Eq. ([34) shows that the Kirkwood relation p| becomes 



the three-particle correlation functions factorize into the 
product of pair correlation functions. This remarkable 
property continues to hold for higher-order correlation 
functions: The same calculation as above shows that 



n 



1 



n fc-2 



fc-1 

II<' 



(35) 



We could have written the denominator on the right-hand 
side of (35) as p k ~ 2 = Il2<a<fc-i( rti -^ so tna * u ^e in 



((1 - n n )(l - n l2 )(l - n l3 )) 
= ((I - n n )(l - n t2 )) ((I - n t2 )(l - n l3 )) 
(1 ~n i2 ) 

This and other similar results are a bit more cumbersome 
than the corresponding expressions in the low-density re- 
gion when one expresses them in terms of particle cor- 
relation functions; in terms of the vacancy correlation 
functions, however, they are as simple. 
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Standard correlation functions do not exploit the na- 
ture of the steady states characterizing the RP. For the 
RP it makes sense to compute the density of islands of 
vacant sites (in the low-density region). The density E s 
of islands of s empty sites delimited at both ends by par- 
ticles can be written as (no(l — n\) ... (1 — n s )n s+ i) in 
terms of the standard correlation functions, yet it is eas- 
ier to compute this quantity directly. Indeed, we need to 
count the total number of configurations of the type 

N—2 particles 



• ooo ooo • ojojojo|ojojo|o (39) 



s vacant li s 



V — s vacancies 



Two particles delimiting the island of s empty sites are 
explicitly presented in (39), possible locations of other 



particles are shown by |. Thus 



E s = 



N 

T 



(V-s-U 
V N-2 J 



IV- 
\N- 



1 - 2p 
1- P \l-P 



(40) 



It is also clear that neighboring islands are uncorrelated, 
e.g. the density E Sl S2 of the two-island pattern 

• ooooooo»oooooooo • 



si vacancies 



S2 vacancies 



is E Su 
is similarly 



p 1 E Sl E S2 . The density of p— island patterns 



E x 



-O-i) 



(41) 



a=l 



DIFFUSION AND SELF-DIFFUSION 



In this section we establish our chief analytical re- 



sult (12) for the density-dependent diffusion coefficient. 



We also derive ( 14 ) which gives the amplitude of self- 
diffusion. Wc consider the symmetric version: Hopping 
to the left and right occur with the same unit rate. 

It is generally impossible to compute D{p) for interact- 
ing lattice gases. The symmetric simple exclusion process 
(SSEP) is a rare exception: In this case the averages 
( n k(t)} satisfy a closed equation which on the hydrody- 
namic scale reduces to Eq. (11) with D = 1. For inter- 



acting lattice gases the governing equations for the aver- 
ages involve higher-order correlation functions, and this 
hierarchical nature prevents the derivation of the diffu- 
sion coefficient. To circumvent this obstacle, one can try 
to compute the diffusion coefficient via the Green-Kubo 
formula [19] . This formula expresses the diffusion coef- 
ficient through the current-current correlation function. 
Current-current correlations are still poorly understood 
for (deterministic) mechanical systems [19] . For (stochas- 
tic) lattice gases these correlations are more tractable [3] . 
The Green-Kubo formula can be schematically written in 



the form 



D(p) 



H(p)- / dtC(t) 



where H (p) is the average total hopping rate and 



CO 

t=-oo 



(42) 



(43) 



is the compressibility (more precisely, the latter is pro- 
portional to x). An explicit expression for C(t) can be 
found in Ref. [3J. This integral contribution has never 
been computed, apart from a few cases where it has been 
proven to be equal to zero [3] . Fortunately, it happens to 
be the case for the RP (see Appendix [Bj), so we can drop 
this term. 

The average total hopping rate from a site in the sym- 
metric version is two times larger than the current in the 
asymmetric version. This gives 



D{p) 



x(p) 



(44) 



which allows us to compute the diffusion coefficient. Us- 
ing ([§ and ([9} we determine the compressibility 



X (p)=p(l-p)\l-2p\ 



(45) 



Plugging ^ and (45) into (44) we arrive at (12). 

Consider now the process of self-diffusion, that is, we 
tag a particle and follow its trajectory. The lattice gas 
is assumed to be at equilibrium and the tagged particle 
is indistinguishable from other particles. Without loss 
of generality, we can assume that the tagged particle is 
initially at the origin. Generally for an arbitrary lattice 
gas with symmetric hopping and in an arbitrary spatial 
dimension d, the average position of the tagged particle 
does not change with time 



<X(i)> = 



(46) 



while the mean-square displacement exhibits a diffusive 
growth 



(X 2 (t)} =2dD s {p,d)t 



(47) 



The coefficient of self-diffusion D s (p, d) is unknown even 
for the SSEP when the spatial dimension is d > 2. One 
should also keep in mind that the diffusion of the tagged 
particle in higher dimensions is generally described by 
the self-diffusion matrix, so one should replace (47) by 
matrix generalization. 

The coefficient of self-diffusion may vanish in one di- 
mension. This happens for the SSEP where the mean- 
square displacement was found [2"0Tf2"5"] to exhibit a re- 
markable sub-diffusive growth, (A 2 ) = Dssep(p)"V% 
with 



SSEP 



(P) 



2 1 



(48) 



As for normal self-diffusion, P(X,t) = Prob(X(i) = X) 
is a Gaussian distribution [24] . The anomalously slow \ft 
growth is caused by the fact that the original particle or- 
dering is forever preserved for the SSEP in one dimension. 
The same growth is valid for other exclusion processes 
with symmetric hopping satisfying two constraints: 

1. No more than one particle per site. 

2. Only nearest-neighbor jumps are allowed. 

The amplitude D(p) generally depends on the de- 
tails of the process. The derivation of (48) given in 
Refs. [221 [5S] suggests that the amplitude ©(p) can be 
expressed through the diffusion coefficient D(p) and the 
compressibility x(p) : 



D(p) = 



2 x(p) 
\fv P 2 



(49) 



Substituting (12) and (45 1 into (49) we arrive at (14) 



Comparing the RP and the SSEP we see that the self- 
diffusion is faster in the latter model, D(p) < Dssep(p); 
in contrast, the relaxation to the equilibrium proceeds 
faster for the RP since D(p) > -Dssep = 1. 

We emphasize that even for the one-dimensional SSEP, 
the normal diffusion behavior (47) holds in finite systems, 



e.g. on the ring or in the case of open boundary condi- 
tions. In these settings, the self-diffusion constant has 
been computed (see [26 for one such calculation) and it 
was found to scale as L^ 1 with system size, so it van- 
ishes in the thermodynamic limit. Note also that the 
two growth laws, \ft and L~ 1 t : match at t c ~ L 2 which 
is the usual diffusion time. The self-diffusion coefficient 
was also computed [2"TH2"9"] for the asymmetric SEP in 
finite systems where it was found to scale as L" 1 / 2 with 
system size, reflecting the (X 2 ) c ~ t 2 l z growth law of the 
variance for infinite systems. 



VI. LATTICE GASES WITH AN EXTRA 
NEXT-NEAREST-NEIGHBOR REPULSION 

The RP studied in the previous sections was based on 
the zero-temperature dynamics compatible with Hamil- 
tonian ([I]) accounting for the nearest-neighbor repulsion. 
In this section we modify the Hamiltonian to include the 
repulsion with next-nearest-neighbors, namely we con- 
sider the RP with Hamiltonian The coupling con- 
stants are assumed to obey J\ > J2 > 0, so that both 
interactions are repulsive and the strength decays with 
separation. The precise values of the coupling constants 
are immaterial, the above assumptions suffice to devise 
a unique zero-temperature dynamics. The general prop- 
erty of zero-temperature dynamics is that an attempted 
hop to an empty site is performed if the energy is not 
raised and rejected otherwise. Thus the dynamical rules 
can be summarized as follows: 



1. If an attempted hop would decrease the number 
of nearest-neighbor pairs of particles, it is always 
performed. 

2. If the number of nearest-neighbor pairs remains the 
same, and the number of next-nearest-neighbors 
would not increase, the hop is always performed. 

3. Otherwise, the attempted hop is never performed. 



A. Steady states and current 

The system wanders on a phase space of maximal- 
island configurations, all of them occurring with the same 
probability. Admissible maximal-island configurations 
are now different in the regions < p < |, 3 < p < |, 
I < p < §, and I < p < 1. It suffices to examine two 
low-density regions: p < | and | < p < |. 

When p < i, the steady states are maximal-island 
configurations with islands of vacant sites of length > 2: 

• oo»ooo»oooo»oo»oooo 

The total number of configurations of type 

N— 1 doublets 
, s 

• 00 |o|o|o|o|o|o|o| o (50) 

S v ' 

V— N — 2 vacancies 

Q, denotes a possible location of the doublet • o) is equal 
to ( n^ 1 )- ^he total number of admissible maximal- 



island configurations is therefore 



e = 



L V - N - 1 



N 



N - 1 



P < 



(51) 



where the -4 factor allows us to consider configurations 
(50) starting with the particle. To compute the current 



(in the asymmetric version) we merely count the number 
of configurations 



JV— 1 doublets 



•000 |o|o|o|o|o|o|o| o 



V— N— 3 vacancies 



and divide by the total number of configurations (51) 
We obtain 



J 



N ( V N N i 2 ) N{V-2N) 



(^r 1 ) Uy-N-i) 



(52) 



When I < p < 5, admissible maximal- island configu- 
rations have islands of vacant sites of length 1 or 2: 



• oo»oo»o»oo»oo»o»o» o 



(53) 



Counting the total number of configurations of type 

V — N vacancies 

, " s 

• o |»o4,»o4.»o4.»o4,»o4,»o|»o| (54) 
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yields { V _ N ) • [In (54]), we denote by J, a possible location 
of a vacancy] The total number of admissible maximal- 
island configurations is therefore 



L ( N 
N\V-N 



\<P<\ 



(55) 



in agreement with the general prediction ( 20 ) . To com- 



pute the current we count the number of configurations 



V— N— 1 vacancies 



• o • o o • o |»o|»o|»o|»o|»o4,»o4. 



and divide by (|55|) to yield 
J 



N (v-n-i) (2N-V)(V-N) 



L 



t N \ 
W-N) 



L(N - 1) 



(56) 



Taking the thermodynamic limit of ( 52 ) and ( 56 1 , ex- 
tending to the high-density region, and collecting all the 
results we arrive at 



J(p) 



K 1 


-3p) 




1- 


-2p 




(3p 


-f)(l- 


2p) 


(2p 


P 

-1)(2- 


3p) 




1-p 




(3p 


-2)(1- 


P) 


2p-l 



< p < ! 

1 < P < 1 



(57) 



This agrees with the announced general result (21) spe- 
cialized to m = 2. 



B. Correlation functions 

In the p < \ region, the pair correlation is easily com- 
puted when the separation does not exceed 2: 

( n i) = P-, ( n i n t+i) = (nin i+2 ) = (58) 

To compute {riirii+i) when t > 3, we first consider 
configurations with q — 1 particles and £ — q vacancies 
between sites i and i + I: 



q—l doublets 



N—q — 1 doublets 



o • oo|ojojojotoojo|ojojojojojoj 



-2t/+l vacancies 



V— N —£-\-2q— 2 vacancies 



Particles at sites i and i+£ are explicitly presented; other 
particles belong to • o doublets, possible locations of dou- 
blets are shown by The total number of such configu- 
rations is 



- 2q - 1\ (V - N - I + 2q - 1 
q-l ){ N-q-1 

Summing over q from q = 1 to q = \_£/3\ we get 

(V-N-i+2q-l 



(59) 



[t/3} 

E 

9=1 



{niUi+t) = ~ Yl V q _ i 



which in the thermodynamic limit becomes 

, , l ^ [£-2q-l\ p^l-Sp)^ 

{nm+e) = E i , - 1 J (i - 2p)^ (60) 

In the J < p < I region analogous computations give 



[1/2] 

E 



(1 - 2p) l - 2 i 
£ - 2qJ p^il-Zpf-^ 



(61) 



For higher-order correlation functions, Eqs. ( 34 1— ( 35 ) 



continue to hold. This is valid in the entire density range. 



C. Compressibility 



When p < 1 , we use the definition ( 43 ) together with 



Eq. ( 58 ) and write the compressibility in the form 



X = p(l - P) - 4p 2 + 2^T [(rum+t) - p 2 ] 



e>3 



To compute the sum in (|62j) we generalize it to 

s (p, A ) = E [( n i n i+e) - p 2 ] xe 



(62) 



(63) 



e>3 



and then deduce the necessary sum by taking the A j" 1 
limit. The advantage of this trick is that we can split the 
sum ( 63 1 into two sums which are both converging when 



A < 1. One gets (see Appendix |C|) 

X (p)=p(l-2p)(l-3p), p<\ (64) 

Similar calculations (see Appendix [C]) give 

X (p)=p(l-2p)(3p-l), \<p<\ (65) 

Using mirror symmetry we then determine x{p) m the 
entire density range: 



fp(l-3p)(l-2p) 
xW = Jp(3,-l)(l-2p) 



< p < | 
3 <P< 5 



< 



I (l-p)(2p-l)(2-3p) ±<p 
[(l-p)(2p-l)(3p-2) |<p<l 



D. Diffusion and self-diffusion 



(66) 



diffusion coefficient 



D{ P ) 



Using (44 1 together with (57 1 and (66) we compute the 



(l-2p)- 2 
p- 2 

a-p)- 2 

(2p-l)- 2 



< p < I 

\<p<\ 

1 2<P<1 

1 <p< 1 



(67) 
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in agreement with the announced result ( 22 ) specialized 
to m = 2. 



We can apply the Green-Kubo formula (44) since the 
gradient condition holds; it does not hold manifestly (as 
for the SSEP), but it holds at equilibrium. This can be 
verified by straightforward analysis. 



The variance of the tagged particle is described by ( 48 1 . 
Combining the general relation (49) with (66) and (67) 
we compute the amplitude in ( 48 ) : 



2>(P) 



ITT yCT 



p(l 


-3p) 




< p < I 


(3p 


-1)(1- 


2p) 


\<P<\ 


(2p 


-1)(2- 


3p) 


3 <P< 3 


(1- 


p)(3p- 


2) 


§<P<1 



(68) 



VII. GENERALIZED REPULSION PROCESSES 

Consider now the generalized repulsion process (GRP) 
characterized by the Hamiltonian (15), with coupling 



constants satisfying the set of constraints (16). One can 
verify the structure ( 18 ) of admissible maximal-island 



configurations in regions ( 17 1 by extending the arguments 



given for the RP. One can also check that configurations 
occur with equal probabilities; these probabilities depend 
on the region and on the density, and are equal to the in- 
verse of the total number of configurations. The deriva- 
tions of the announced results for the total number of 



configurations, Eqs. ( 19 )— ( 20 ), and the steady state cur- 
rent |2ll) are also straightforward. The expression (22) 



for the diffusion coefficient then follows from the Green- 



Kubo formula ( 44 ) when we combine it with expressions 



(21 ) for the current and with formula 



p[l - (to + l)p] [1 - mp] 
p[(k + l)p-l][l-kp] 



< p < 



i 

m+l 
<P<\: 



(69) 



fe+i 



giving the compressibility; k — 2, . . . , m in Eq. (69). To 
justify the applicability of the Green-Kubo formula (44) 
we use again the gradient condition; checking of the va- 
lidity of the gradient condition at equilibrium is tedious, 
but straightforward. 

To complete the derivation we must establish Eq. ( 69 1 . 



The approach that has been used previously is still ap- 
plicable: One determines the connected pair correlation 
function and then one finds the compressibility via the 
definition (43). A long calculation was necessary even 
when to = 2. To circumvent the necessity for such a 
calculation, we employ a different strategy. The com- 
pressibility can be found (see e.g. [B]) from relation 



d 2 F 

Here F is a free energy per lattice site: 



F 



lim — In (2 

L—t-oc L 



(70) 



(71) 



The computation of the total number of configurations 6 
is much simpler than the computation of the connected 



pair correlation function, and then from Eqs. ( 70 )— ( 71 ) 
we directly deduce the compressibility. 
For the RP, we get 

F = plnp - (1 - p) ln(l - p) + (1 - 2p) ln(l - 2p) 

in the low-temperature region. Using this result in con- 
junction with (70) we re-derive (45) thereby providing a 
useful check of the consistency of our results [Eqs. ^ 
and ^] for the connected pair correlation function. 

Similarly for the RP involving next-nearest-neighbor 
repulsive interactions, we use (|5l|) and (55) to compute 



F = p In p - (1 - 2p) ln(l - 2p) + (1 - 3p) ln(l - 3p) 
for < p < | and 
F= (1- 2p) ln(l - 2p) -plnp + (3p - l)ln(3p- 1) 

Using these results in conjunction with 



for i < p < 



( 70 1 we re-derive ( 66 ) 



Generally, one can use Eqs. ( 19 )— ( 20 ) to establish the 



free energy. This leads to Eq. ( 69 i for arbitrary to 



VIII. DISCUSSION 

We studied dynamical behaviors of one-dimensional 
stochastic lattice gases. We assumed that each site is 
occupied by at most one particle and that particles in- 
teract through repulsive forces whose span can be arbi- 
trary. We endowed the system with a zero-temperature 
dynamics, so that the hops which would have led to the 
increase of energy are forbidden. When the strength 
of interactions monotonically and rapidly decreases with 
the separation between the particles, namely the cou- 
pling constants satisfy the set of constraints (16), inter- 
actions can be treated in a lexicographic order. In the 
case of symmetric nearest-neighbor hopping, we analyt- 
ically determined the density-dependent diffusion coeffi- 
cient. We also computed the variance of the displacement 
of a tagged particle. 

The repulsion process (RP) with nearest-neighbor in- 
teractions is related to a few other lattice gases. One such 
model is the facilitated exclusion process [3U] which is 
well-defined in the high-density regime and has the same 
structure of steady states and the same current as the 
RP, even though away from the steady state the dynam- 
ical rules of the two processes differ. Another is the KLS 
model [10 which reduces to the RP at a certain point of 
the parameter space. Our the results for the current can 
be extracted from the earlier work [T^] on the KLS model, 
and they apparently were known much earlier [3"Tj . The 
correlation functions, the density-dependent diffusion co- 
efficient, and the amplitude characterizing self-diffusion 
haven't been investigated even in the context of the RP 
with nearest-neighbor interactions. Most of our results 
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for the RP with nearest-neighbor interactions apparently 
admit a generalization to the KLS model. The tools 
are different, namely one employs a transfer-matrix tech- 
nique [HI H3] ; the emerging results are very similar [T5] . 

The vast majority of lattice gases which have ever 
been studied involve only nearest-neighbor interactions. 
Therefore the generalized repulsion processes (GRP) 
with interactions whose span can be arbitrary form an 
interesting class of lattice gases with long-ranged inter- 
actions which, remarkably, admit the exact analyses. 

In this work we limited ourselves to the GRP on the 
ring. It would be interesting to study the GRP with 
asymmetric hopping in an open setting. Interesting be- 
haviors have been reported for lattice gases with the cur- 
rent J(p) having one or two maxima 1 1 2] . For the 
GRP we have potentially many maxima, viz. 2m max- 
ima for the system with Hamiltonian (15), which might 
result in a rich phase diagram. 

Another avenue for future work is to look for large devi- 
ations [32] . Interesting singular behaviors of the large 
deviation function have been recently found }14| in the 
context of the KLS model, in the point in the parameter 
space close to the RP. Similar singular behaviors proba- 
bly arise for the RP, and perhaps more subtle behaviors 
will appear for the GRPs. In Appendix [D] we specialize 
the macroscopic fluctuation theory |32j describing large 
deviations in diffusive lattice gases to repulsion processes. 

Finally we note that it would be interesting to perform 
an analytical investigation of the driven GRP in the set- 
ting with open boundaries. The influence of repulsion can 
be striking, e.g. very rich phase diagrams were found in 
similar models (see [T2J |33j and references terein) . 



Acknowledgments 

I am grateful to Uttam Bhat, Kirone Mallick, Baruch 
Meerson, and Sid Redner for discussions. 



Appendix A: Higher-order connected correlation 
functions 



Here we establish explicit expressions for the connected 
correlation functions (12345) c and (123456) c . We then 
compute the number of exponential terms arising in the 
p— particle connected correlation function for arbitrary p. 

Let's begin with already known connected correlation 
functions given by Eqs. pi), (37), and (38). It is useful 



to think about terms in these equations as representing 
contributions from the diagrams: 12 in the case of pi); 



123 in the case of (37 1; 12 34 and 1234 in the case of (38 1 



Continuing this line of reasoning we see that diagrams 



contribute to (12345) c = (n^n^n^n^n^) c . Therefore 

(12345) c = AR l2 -' ll+i '-- i3 + AR 1 "- 11 ^-^ 

BR iB ~ il , R = -^— (Al) 



1-P 

The amplitudes in front of the terms in the top line of 



Eq. (Al) are equal each other due to symmetry. The 



ratio R = lies in the range — 1 < R < in the 

l— p ° 

low-density region p < | . 
Similarly, the diagrams 

12 34 56 123 456 123456 1234 56 123456 
contribute to (123456) c . This leads to 

(123456) c = f7^-i 5 +i4-i 3 +i2-ii + y R i 6 -U+i 3 -ii 

+ WR l6 ~ 11 (A2) 



The amplitudes appearing in the second line of Eq. ( A2 ) 



are equal due to symmetry, and they coincide with one 
of the amplitudes in the first line (we don't know a con- 
ceptual explanation for this equality). 



To determine the amplitudes A and B in Eq. ( Al ) we 
set 1 = 2 (that is, i\ — 12). Using the definition (35) 
together with relation 



(n-pf 

we find 

(22345) c 



n 2 — 2np + p 2 = n — 2np + p 2 
(l-2p)(n-p)+p(l-p) 



(A3) 



= (l-2p)<2345) c + p(l-p)<345) t 

= (l-2p)[p(l- p )fR^-i2+iB-i4 

+ (1 - 2pfp{\ - p)R 1 ^ 

+ (l-2p)[p(l-p)] 2 R^ 



On the other hand, Eq. ( Al ) reduces to 

(22345) c = AR i5 - i3 + AR l3 - l2+is - i4 + BR 5 ~ i2 
when i\ = ii. Comparing these two results we get 
A = (1 - 2p)[p{\ - p)} 2 , B = (1 - 2p) 3 p(l - p) 



Similarly we set 1 = 2 in Eq. ( A2 ) and reduce it to 
(223456) c = UR i6 ~ i5+iA ~ ia +VR ie ~ iA+i3 ~ i2 

+ WR i6 ~ 12 (A4) 



12 345 



123 45 



12345 



On the other hand, using |35|) and (A3) we obtain 
(223456) c = (l-2p)(23456) c + j o(l- j o)(3456) c 

= (1 - 2p)A (iJ*6-<4+*3-i2 + R i6-h+U-i2^ 
+ (1 - 2p)BR 6 - 12 + [p{\ - /9 )] 3 ii i 8-i5+i4-i3 

+ (l-2p) 2 [p(l-p)] 2 R^ 
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Comparing this with (A4) we determine the amplitudes 



U=[p{l-p)f 
V=(l-2p) 2 [p(l-p)] 2 
W=(l- 2pfp{\ - p) 

Consider now the general case and denote by D p the 
number of diagrams (equivalently, the number of expo- 
nential terms) appearing in the p— particle connected cor- 
relation function. The number of diagrams starting with 
12 is equal to D p ^ 2 , the number of diagrams starting 
with 123 is equal to D p _ 3 , etc. Therefore 

D p = D p -2 + D p -3 + . . . + D 2 + 1 

which can be re-written as the recurrence 



D p = Dp_ 2 + L>p_ 



(A5) 



defining the Fibonacci numbers. The standard initial 
condition for the Fibonacci numbers is F\ — F 2 = 1, 
while in our case D 2 = D3 = 1. Hence D p = F p -i. 



Appendix B: Repulsion processes and gradient 
condition 



Zero-temperature dynamics which are compatible with 
Hamiltonian ([!]) and involve only nearest-neighbor hop- 
ping proceed according to the following rules: 

1. If a hop (to a neighboring empty site) would result 
in decrease of the number of nearest-neighbor pairs 
of particles, AH = — 1, it is always performed. 

2. If a hop would cause neither decrease nor increase 
of the number of nearest-neighbor pairs of particles, 
AH = 0, it is performed with probability p. 

3. If a hop would increase of the number of nearest- 
neighbor pairs of particles, AH = 1, it is never 
performed. 

The probability p is a parameter of the model. Three 
values of p play a special role: 

• p = 0. Only strictly energy lowering hops are al- 
lowed. In this case, the system quickly gets trapped 
in a jammed configuration. 

• p = \. In zero-temperature spin-flip dynamics, 
this special value is prescribed by Glauber dynam- 
ics, and it is known to simplify analysis in a few 
tractable cases [9]. 

• p = 1. This value is most efficient in simulations, 
and it is typically used in Metropolis and heat-bath 
algorithms. 



As long as p > 0, its value usually plays a minor role. 
In this paper we have chosen p = 1 although, similarly 
to zero-temperature spin-flip dynamics, the choice P = \ 
is theoretically advantageous. To see why this is so, let 
us consider the current Ji^+i through the bond (i, i + 
This current is fully determined by occupation variables 
Tii—x, Hi, n i+ i, n i+2 . When p = ~, the current reads (we 
assume that the hopping is symmetric): 



2Ji 



M+i 



H+l 



+ (tli-i 



H+2 



H+l) 



(Bl) 



while for p — 1 one gets 



Ji,i+l = rij_inj(l - n i+ i) - (1 - ni)n i+1 n i+2 

+ (m - n. i+ i)(l - ni_i)(l - n i+2 ) (B2) 

The current corresponding to the P — \ case is not 
merely expressed by a more compact formula than the 
current corresponding to p = 1. A crucial property is 
that for p — I the current has a gradient form. Indeed, 
massaging (Bl ) we re- write it as 



2Jm+i = N ~ n i+i\ + [rii-in i+ i - n t n i+2 ] 

+ [rii-iUi - n t n i+1 ] + [mn i+1 - n l+1 n l+2 ] 
- 2[n t ^in l n t+ i - nin i+ xni+%\ 

where each term in the square brackets is a (discrete) 
gradient of some local function of n = {rij}. For p = 1, 
and for any p ^ | , the current does not have a gradient 
form. For instance, for p = lwe can re-write (B2| as 



Ji,i+l = [th - n l+ i] + [ni-iUi+i - riin l+2 ] 
- [ni-in z n i+ i - nini+mi+2} 
+ ni-\nin i+ 2 - ni-in i+ in i+2 

Each term in the square brackets is a gradient; the terms 
in the last line cannot be converted into a gradient form. 

If the current across the bond admits a gradient rep- 
resentation, the computation of the diffusion coefficient 
greatly simplifies, namely the integral term in Eq. (42) 
vanishes [3], and we can a simple formula (44 1. This 



formally justifies the results of Sect. [Vjfor the RP with 
p = i. We now recall that the hydrodynamic regime is 
very close to the (local) equilibrium where energy low- 
ering hops are no longer possible. Therefore there is no 
distinction between the dynamics with any p > 0, all hops 
occur with the same rate p. Therefore the expression for 
the diffusion coefficient is actually the same (apart from 
the trivial rescaling, D — > pD). 



Appendix C: Compressibility 



Let us split the sum ( 63 ) into two sums 



S(p, A) = J2{n inm ) X e -p 2 J2^ = Si- S 2 



f>:s 



e>3 
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To compute Si we use (60), write t = 3q + m, change the 
itiorj 

pX^_ 

l-'3 P ]- q 



order of summation, and after a bit of massaging we get 



Si 



p 



9>1 v r/ m>0 

P x 3 ^ q - 



l-3p 



1 - 2p 



g>l v 



2p 



1 — A 



1 - 2p 



1 - 2p - A(l - 3p) - pA 3 
Therefore 

p 2 A 3 



5(p,A) 



2 A 3 



1 - 2p - A(l - 3p) - pA 3 1 - A 



from which we deduce liniA->i iS(/9, A) = 3p 3 . Plugging 



this into (62) we arrive at the announced result (64). 

Similarly one can determine x m the | < p < \ region. 
As previously, it suffices to compute 



p(l- P )+2S 1 (p,X) 



2p 2 X 
l-X 



(CI) 



where Si — ^2i >2 { n i n i+i) A • We plug (61 1 into Si, write 
£ = 2q + m, and change the order of summation to yield 



^ pi- 1 ^ \m 

q>l r m=0 v 

= J2^ x ( i - 2 p)+ 3 p- i y 

q>l P 



-P- 



[A(l-2p)r(3p-l)«- 



p-A 2 [A(l-2p) + 3p-l] 



Substituting this result into (|C1| and taking the A j" 1 
limit we arrive at the announced result ( 65 ) . 



Appendix D: Large deviations 

The macroscopic fluctuation theory (MFT) has been 
recently developed to probe large deviations in lattice 
gases [32]. Mathematically, one must solve two coupled 
partial differential equations 



dtp = d x [D(p) d x p] - d x [a(p) d x p] 
d t p = -D(p)d xx p - \& '{p){d x pf 



(Dl) 



for the density field p(x, t) and the conjugate momentum 
field p(x,t). These equations contain the diffusion coef- 
ficient D(p) and the quantity a(p) which characterizes 
the variance of the flux. More precisely, consider a one- 
dimensional system of a very large length L with open 
boundaries which are in contact with reservoirs of parti- 
cles with density p. The variance of the the integrated 
flux gives <r(p): 



t— »oo t 



L 



(D2) 

The equilibrium origin of the quantities D(p) and <r{p) is 
emphasized by relation 



d 2 F(p) 2D(p) 



dp 2 



(D3) 



which follows 3, 8 from the fluctuation-dissipation the- 
orem. Combining (D3| with (70) we arrive at a = 2D\, 
which simplifies to 



a{p) = 2J(p) 



(D4) 



for all GRP. Interestingly, a(p) describes the symmet- 
ric version of the process, while J(p) is the steady state 
current in the totally asymmetric version. We have com- 
puted J(p), see ^ and (21 ), and hence we know a(p) for 
an infinite class of repulsion processes. 



[1] R. J. Baxter, Exactly Solved Models in Statistical Me- 
chanics (Academic Press, London, 1982). 

[2] G. M. Bell and D. A. Lavis, Statistical Mechanics of Lat- 
tice Models (Wiley, London, 1989). 

[3] H. Spohn, Large Scale Dynamics of Interacting Particles 
(New York: Springer- Verlag, 1991). 

[4] C. Kipnis and C. Landim, Scaling Limits of Interacting 
Particle Systems (Springer, New York, 1999). 

[5] B. Derrida, Phys. Repts. 301, 65 (1998). 

[6] G. M. Schiitz, Exactly Solvable Models for Many-Body 
Systems Far From Equilibrium, in Phase Transitions and 
Critical Phenomena, Vol. 19, eds. C. Domb and J. L. 
Lebowitz (Academic Press, London, 2000). 

[7] R. A. Blythe and M. R. Evans, J. Phys. A 40, R333 
(2007). 

[8] B. Derrida, J. Stat. Mech. P07023 (2007). 
[9] P. L. Krapivsky, S. Redner, and E. Ben-Nairn, A Kinetic 
View of Statistical Physics (Cambridge University Press, 



Cambridge, 2010). 
[10] S. Katz, J. L. Lebowitz, and H. Spohn, J. Stat. Phys. 34, 
497 (1984). 

[11] J. Krug, Phys. Rev. Lett. 67, 1882 (1991). 

[12] J. S. Hager, J. Krug, V. Popkov, and G. M. Schiitz, Phys. 

Rev. E 63, 056110 (2001). 
[13] R. K. P. Zia, J. Stat. Phys. 138, 20 (2010). 
[14] G. Bu nin, Y. Kafri, an d D. Podolsky, EPL 99, 20002 

(2012); |arXiv:1301. 17081 
[15] P. L. Krapivsky, in preparation. 
[16] P. L. Krapivsky and J. Olejarz, [arXiv:1303.5128"l 
[17] J. G. Kirkwood, J. Chem. Phys' 3, 300 (1935). 
[18] R. Balescu, Equilibrium and Nonequilibrium Statistical 

Mechanics (Wiley, New York, 1975). 
[19] P. Resibois and M. De Leener, Classical Kinetic Theory 

of Fluids (Wiley, New York, 1977). 
[20] T. E. Harris, J. Appl. Prob. 2, 323 (1965). 
[21] D. G. Levitt, Phys. Rev. A 8, 3050 (1973). 



14 



[22] P. M. Richards, Phys. Rev. B 16, 1393 (1977). 
[23] S. Alexander and P. Pincus, Phys. Rev. B 18, 2011 
(1978). 

[24] R. Arratia, Ann. Probab. 11, 362 (1983). 

[25] H. van Beijeren, K. W. Kehr and R. Kutner, Phys. Rev. 

B 28, 5711 (1983). 
[26] T. Antal, P. L. Krapivsky, and K. Mallick, J. Stat. Mech. 

P08027 (2007). 

[27] B. Derrida, M. R. Evans, and D. Mukamel, J. Phys. A 

26, 4911 (1993). 
[28] B. Derrida, M. R. Evans, and K. Mallick, J. Stat. Phys. 

79, 833 (1995). 



[29] B. Derrida and K. Mallick, J. Phys. A 30, 833 (1997). 
[30] A. Gabel, P. L. Krapivsky, and S. Redner, Phys. Rev. 

Lett. 105, 210603 (2010). 
[31] H. Brandstetter, Diploma thesis, Univ. Munich, 1991 

(unpublished). 

[32] L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, and 
C. Landim, Phys. Rev. Lett. 87, 040601 (2001); ibid 94, 
030601 (2005); J. Stat. Mech. (2007) P07014. 

[33] M. Dierl, P. Maass, and M. Einax, Phys. Rev. Lett. 
108, 060603 (2012); M. Dierl, M. Einax, and P. Maass, 
larXiv:1304.5953l 



